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1. 


This  final  report  is  being  written  to  briefly  describe 
the  research  conducted  under  contract  No.  F19628-79-C-0027 
from  December  1,  1978  to  October  31,  1981.  The  research 
has  been  carried  ou  the  theme  of  the  contract  title:  The 
Study  of  Factors  Affecting  the  Precise  Estimation  of  the 
Earth's  Gravity  Field.  The  research  has  been  conducted  by 
three  graduate  research  associates,  three  post  doctoral 
researcher,  and  the  principal  investigator.  The  research 
results  have  been  reported  in  detail,  in  12  Scientific  Reports 
of  the  contract  and  8  papers  presented  at  various  meetings. 

In  addition,  quarterly  reports  have  been  made  to  the  Contract 
Manager  from  the  start  of  the  project. 

It  does  not  seem  appropriate  to  detail  the  research 
that  has  been  described  earlier.  However,  it  is  appropriate 
to  describe  areas  of  study,  the  reports  associated  with  each 
one,  and  the  most  significant  contributions. 

A  complete  list  of  reports,  by  scientific  report  number 
is  given  in  section  4.  A  list  of  the  papers  presented  at 
meetings  is  given  in  section  5,  and  a  list  of  persons  employed 
under  the  contract  is  given  in  section  3. 


Areas  of  Studj 


The  research  can  be  described  in  the  following  four 
categories:  new  analysis  techniques  for  data  given  on  a 

spherical  grid;  analysis  of  satellite-to-satellite  tracking 
data;  global  gravity  field  modelling;  and  the  vertical 
datum  problem. 

New  Analysis  Techniques  For  Data  Given  on  a  Spherical  Grid 


Many  geodetic  problems  require  the  analysis  of  data 
given  on  a  spherical  surface  representing  the  earth.  As 
an  example,  we  have  the  calculation  of  geoid  undulations 
and  potential  coefficients  from  mean  gravity  anomalies.  As 
the  number  of  surface  elements  increases,  and  the  number 
of  quantities  of  interest  increase  the  computational  effort 
can  significantly  increase.  This  is  especially  true  if  an 
optimal  estimation  of  the  parameters,  considering  data  noise 
and  other  factors,  is  to  be  done. 

Colombo  (Sci.  Rpt.  No.  1,  1979)  examined  the  application 
of  least  squares  collocation  to  data  given  on  a  uniform 
spherical  grid.  In  this  report  a  theoretical  analysis  was 
carried  out  with  both  point  and  area  mean  data.  Of  special 
theoretical  interest  was  the  spherical  harmonic  analyses 
of  the  data.  The  techniques  discussed  in  this  report  were 


more  thoroughly  analyzed  by  Colombo  (Sci.  Rpt.  No.  7,  1981). 

In  this  report  the  role  of  two-dimensional  Fourier  analysis 
in  spherical  harmonic  analysis  is  discussed.  These  procedures 
were  used  to  develop  highly  efficient  computer  programs  to 
generate  data  on  a  grid  from  a  given  set  of  spherical  harmonic 
coefficients,  and  to  generate  coefficients  from  a  set  of 
given  data,  either  point  or  mean  values.  As  an  example, 
the  computer  time  (on  an  Amdahl  470/V6,  II)  to  generate  a 
set  of  spherical  harmonic  expansion  to  degree  300  from  64800 
l°x  1°  mean  values  was  36  seconds. 

This  Colombo  report  was  also  important  because  it  showed 
how  the  size  of  the  mean  blocks  being  used  as  data,  aliases 
the  coefficients  being  computed.  This  was  crucial  in  identi¬ 
fying  sampling  error  as  well  as  data  noise  in  estimating 
the  accuracy  of  the  computed  quantities. 


Analysis  of  Satellite-to-Satell ite  Tracking  Data 


The  main  concern  in  this  area  is  the  development  of 
techniques  for  the  processing  of  SST  tracking  data  that  might 
be  expected  from  the  GRAVSAT  mission.  Rapp  and  Ha.iela  (Sci. 
Rpt.  No.  3,  1979),  first  studied  the  accuracy  to  be  expected 
for  l°x  1°  mean  anomaly  determination  from  a  GRAVSAT  type 
of  mission.  Rummel  (Sci.  Rpt.  No.  6,  1980)  postulated 
various  data  distributions  and  used  least  squares  collo¬ 
cation  to  estimate  the  accuracies  of  various  quantities  of 
interest.  Hajela  (Sci.  Rpt.  No.  9,  1981)  made  the  first 
attempt  at  a  simulation  study  to  examine  the  accuracy  of 
anomaly  recovery  using  collocation.  In  this  study  range- 
rate  data  was  generated  by  Wong  and  Sjogren  given  a  set  of 
l°x  1°  mean  anomalies.  Hajela's  attempt  at  the  recovery  of 
the  starting  anomalies  was  moderately  successful. 

The  fundamental  problem  in  the  above  analyses  was  the 
limitation  imposed  by  data  density  and  the  number  of  data 
points  that  could  be  used  in  a  usual  least  squars  collocation 
solution.  Colombo  (Sci.  Rpt.  No.  10,  1981)  developed  a  collo¬ 
cation  procedure  to  analyze  the  SST  data  to  recover  spherical 
harmonics.  An  actual  simulation  was  run  with  a  6  month  set 
of  data  points  to  determine  the  accuracy  in  which  the  harmonics 
could  be  determined.  This  procedure,  expanded,  perhaps, 
to  consider  orbit  error,  may  be  the  way  for  an  actual  data 
analysis  since  the  method  can  handle  large  data  sets  in  a 
reasonable  amount  of  time  in  an  optimum  way. 

Another  look  at  large  dense  data  sets  has  been  done 
by  Rummel  (1981).  Although  the  method  of  Rummel  was  developed 
for  the  analysis  of  altimeter  data,  its  applicability  can 
be  expanded  to  SST  data. 
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Global  Gravity  Field  Modelling 


We  continued  to  study  ways  in  which  the  earth's  gravity 
field  could  be  determined  and  analyzed  from  existing  data, 

Rapp  (Sci.  Rpt .  No.  2,  1979)  examined  the  problem  of  having 
a  good  anomaly  degree  variance  model  for  covariance  computations. 
New, one  and  two  component  models  were  developed  using  the 
results  of  a  spherical  harmonic  expansion  to  degree  180  made 
in  1978. 

Many  combination  solutions  had  been  carried  out  using 
5°  mean  anomalies.  A  study  was  made  to  see  the  effect  of 
using  5°  or  l°x  1°  data  in  a  combination  solution.  The  results 
were  described  by  Rapp  (Sci.  Rpt.  No.  5,  1980),  who  showed 
that  differences  (between  using  5°  and  1°  data)  of  50%  could 
occur  at  degree  36.  Such  differences  were  later  verified 
in  an  analytical  and  numerical  manner  by  Colombo  (Sci.  Rpt. 

No.  7,  1981). 

In  1981  a  substantial  effort  was  made  to  develop  a  new 
spherical  harmonic  expansion  of  the  earth's  gravity  field. 

This  solution  (Rapp,  Sci.  Rpt.  No.  12,  1981)  was  carried 
out  to  degree  180  using  Seasat  altimeter  data,  new  terrestrial 
l°x  1°  data,  and  a  new  a  priori  set  of  potential  coefficients 
obtained  from  satellite  analysis.  The  combination  procedures 
used  techniques  and  some  programs  developed  by  Colombo  (Sci. 

Rpt.  No.  7,  1981)  earlier. 

Although  our  analysis  has  been  with  spherical  harmonics, 
there  has  been  a  concern  on  the  evaluation  of  functions  at 
the  surface  of  the  earth  from  such  expansions.  Jekeli  (Sci. 

Rpt.  No.  11,  1981)  analyzed  the  downward  continuation  using  some 
very  high  degree  (up  to  degree  300)  expansions.  In  doing  this 
actual  models  of  the  earth's  topography  were  used.  The  results 
obtained  indicated  corrections  to  the  anomalies  and  undulations 
could  reach  86  pgal  and  410  um.  These  values  are  considered 
negligable  with  respect  to  other  error  sources. 

Another  aspect  of  gravity  field  modelling  is  related 
to  the  computation  of  the  gravity  vector  in  space.  Of  special 
concern  is  the  vector,  near  the  surface  of  the  earth,  where 
the  classical  techniques  are  no  longer  valid.  Katsambalos 
(Sci.  Rpt.  No.  8)  studied  this  problem  using  topographic 
models  and  anomaly  fields  generated  by  point  masses.  He 
found  that  classical  techniques  could  be  seriously  in  error 
for  points  on  the  order  of  1  to  10  km  above  the  terrain. 
Katsambalos  developed  significantly  improved  techniques  for 
this  type  of  computation. 
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Vertical  Datum  Problem 


Because  of  sea  surface  topography  the  mean  ocean  surface 
is  not  the  geoid.  Since  height  datums  are  based  on  mean 
sea  level,  such  datums  are  not  defined  with  respect  to  the 
same  equipotential  surface.  The  question  then  arises  as 
to  how  one  generates  a  unique  vertical  datum,  or  corrects 
existing  vertical  datums.  Our  first  study  was  carried  out 
by  Colombo  (Sci.  Rpt.  No.  4,  1980).  In  this  report,  Colombo 
showed  how  it  would  be  possible  to  correct  datums.  The  data 
needs  were  extensive,  requiring  precise  (±10  cm)  station 
coordinates.  Accurate  levelling  between  selected  stations 
and  a  master  station  in  a  network,  and  geoid  (spherical  harmonic) 
and  local  (anomaly)  gravity  data.  With  some  reasonable  assump¬ 
tions  Colombo  showed  the  correction  of  datums  was  a  feasible 
procedure. 

Some  studies  were  attempted  to  determine  the  accuracy 
of  which  a  vertical  datum  correction  could  be  made  with  data 
available  today.  Unfortunately  these  studies  were  inconclusive 
because  of  the  questions  that  exist  on  the  accuracy  and  distri¬ 
bution  of  the  needed  data. 

Rapp  (1980)  also  discussed  the  concept  of  the  vertical 
datum  question.  Here  the  incorporation  of  sea  surface  topography 
and  satellite  altimetry  was  discussed. 
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